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Random skew plane partitions and the Pearcey 

process 

Andrei Okounkov *and Nicolai Reshetikhinl 


Abstract 

We study random skew 3D partitions weighted by q™* and, specifically, 
the q —> 1 asymptotics of local correlations near various points of the limit 
shape. We obtain sine-kernel asymptotics for correlations in the bulk of 
the disordered region, Airy kernel asymptotics near a general point of the 
frozen boundary, and a Pearcey kernel asymptotics near a cusp of the 
frozen boundary. 


1 Introduction 


A plane partition tt = (TTij) is an array of nonnegative numbers indexed by 
(i,j) € that is monotone, that is, 

TTij ^ j_|_s , r, S ^ 0 


and finite in the sense that tt^ = 0 when i + j 3> 0. Plane partitions have 
a obvious generalization which we call skew plane partitions. A skew plane 
partition is again a monotone array (tt^) which is now indexed by points (i,j) 
of a skew shape A//i, where ^ C A is a pair of ordinary partitions. We call p, 
and A the inner and outer shape of tt, respectively. In fact, in this paper we 
will only consider the case when the outer shape A is a a x 6 rectangle. Here is 
an example with p, = (1,1) and a = b = 5. 


7 4 2 1 

5 3 2 0 
6 3 110 
4 110 0 
2 0 0 0 0 


( 1 ) 


Placing TTij cubes over the (z,j) square in X/p gives a three-dimensional object 
which we will call a skew 3D partition and denote by the same letter tt. Its 
volume is |7r| = X] TTy . For tt as in (^1, it shown in Figure Q] 
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Figure 1: A skew 3D partition 


Given a parameter 0 < g < 1, define a probability measure on the set of all 
skew plane partitions with given inner and outer shapes by setting 

Prob(7r) oc g'’"' . (2) 

The corresponding random skew 3D partition model has a natural random 
growth interpretation, the parameter g being the fugacity. Also, a simple bi- 
jection, which should be clear from Figure Q and is recalled below, relates this 
model to a random tiling problem. 

We are interested in the thermodynamic limit in which g —> 1 and both inner 
and outer shapes are rescaled by l/r where 

r = — In g ^ +0 . 

The results of IS] imply the following form of the law of large numbers: scaled 
by r in all directions, the surface of our random skew 3D partition converges to 
a nonrandom surface — the limit shape. This limit shape will be easy to see in 
the exact formulas discussed below. A simulation showing the formation of the 
limit shape is presented in Figure [21 

An important qualitative feature of limit shape is the presence of both or¬ 
dered and disordered regions, separated by the frozen boundary. Furthermore, 
the frozen boundary has various special points, namely, it has cusps (there is 
one forming in Figure (21 it can be seen more clearly in Figures [TT)1 and ITTIl and 
also turning points where the limit shape is not smooth. In Figure [T^HTTl the 
turning points are the points of tangency to any of the lines in the same figure. 

One expects that the microscopic properties of the random surface, in par¬ 
ticular, the correlation functions of local operators, are universal in the sense 
that they are determined by the macroscopic behavior of the limit shape at that 
point. More specifically, one expects that: 
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Figure 2: A large random skew 3D partition 


1. in the bulk of the disordered region, the correlation are given by the in¬ 
complete beta kernel [12 US] with the parameters determined by slope of 
the limit shape (a special case of this is the discrete sine kernel); 

2. at a general point of the frozen boundary, suitably scaled, the correlation 
are given by the extended Airy kernel 11213112; 

3. at a cusp of the frozen boundary, correlation, suitably scaled, are given by 
the extended Pearcey kernel, discussed below and in E3- 

In this paper we prove all these statements for the model at hand. The required 
techniques were developed in our paper Hi, of which this one is a continuation. 
Namely, as will be reviewed below, our random skew plane partition model is 
a special case of Schur process. This yields an exact contour integral formulas 
for correlation functions. The asymptotics is then extracted by a direct albeit 
laborious saddle point analysis. The striking resemblance of the above list to 
classification of singularities is not accidental for, as we will see, these three 
situation correspond precisely to the saddle point being a simple, double, or 
triple critical point. 

We will also see that the frozen boundary is essentially an algebraic curve 
and that it has precisely one cusp per each exterior corner of the inner shape /i. 

We expect that near a turning point the correlations behave like eigenvalues 
of a A: X A: corner of a GUE random N x N matrix, where A^ ^ 0 and k plays 
the role of time. We hope to return to this question in a future paper. 


3 






The results presented here were obtained in 2002-03 and were reported by 
us at several conferences. The period between then and now saw many further 
developments in the field. Most notably, the Pearcey process, which we found 
describes the behavior near a cusp of the frozen boundary arose in the random 
matrix context in the work of Tracy and Widom EDI. Pearcey asymptotics 
for equal time correlations of eigenvalues were obtained earlier by Brezin and 
Hikami EllH and also by Aptekarev, Bleher, and Kuijlaars El- We enjoyed and 
benefited from the correspondence with C. Tracy on subject. 

In El) Ferrari and Spohn derived from the exact formulas of EDI the Airy 
process asymptotics in the case of unrestricted 3D partitions (the /r = 0, a = 
b = oo case in our notation). In a related but technically more involved context, 
the Airy process asymptotics was found by K. Johansson in EDI. 

In EH) the partition function of the random surface model studied here was 
related to the topological vertex of ^ and, thus, to the Gromov-Witten theory 
of toric Calabi-Yau threefolds. They were many subsequent developments, some 
of which are reviewed in m The papers El EH may the closest to the material 
presented here. Also, much more general results on algebraicity of the frozen 
boundary are now available m We would like to thank R. Kenyon and C. Vafa 
for numerous discussions. 

N. R. is grateful to Laboratoire de Physique Theorique at Saclay for the 
hospitality, where part of this work was done and to J.-B. Zuber and Ph. Di 
Francesco for interesting discussions. His work was supported by the NSF grant 
DMS-0070931 and by the Humboldt Foundation. The work of A. O. was par¬ 
tially supported by the Packard Foundation. 

2 Preliminaries 

2.1 Skew 3D partition as a sequence of its slices 

We associate to tt the sequence {A(t)} of its diagonal slices, that is, the sequence 
of partitions 

X{t) = i > max(0,-t), t e Z. (3) 

Throughout this paper, we assume that the outer shape of our skew partition 
is an a X 6 box and, in particular, we will use the letter A to denote diagonal 
slices, not the outer shape. 

Notation A as usual means that A and v interlace, that is, 

Al > /Ti > A 2 > /i 2 ^ A 3 > ... . 

It is easy to see that the sequence {A(t)} corresponds to a skew plane partition 
if and only if it satishes the following conditions: 

• if the slice A(to) is passing through an inner corner of the skew plane 
partition then 

• • • -( A(to ~ 2) ^ A(to ~ 1) ^ A(to) A(to + 1) ^ A(to -I- 2) ; (4) 
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• if the slice A(to) is passing through an outer corner of the skew plane 
partition then 

■■■)>- A(to — 2) A(to — 1) ^ A(to) ^ A(to “f 1) ^ A(to + 2) ^ . . (5) 

For example, the configuration {A(t)} corresponding to the partition is 

(2) ^ (4) ^ (6, 1 ) ^ (3, 1 ) ^ (5, 1 ) ^ (7,3, 1 ) ^ (4,2) ^ (2) ^ ( 1 ). 

We will denote the sequence of inner and outer corners of the inner shape by 
{i^z}i<i<Ar and {Mi}i<i<Ar-i, respectively. We assume that they are numbered 
so that 

Vi < Ui < V2 < U2 < ■ ■ ■ UN-1 < vn ■ 

We also assume that the point f = 0 is chosen so that 

Vi = ■ ( 6 ) 

l<i<N l<i<N-l 


2.2 Connection to tilings 

There is a well-known mapping of 3D diagrams to tilings of the plane by rhombi. 
Namely, the tiles are the images of faces of the 3D diagram under the projection 

{x, y, z) ^ {t, h) = {y - X, z - (x + y)/2). (7) 

This mapping is a bijection between 3D diagrams and tilings with appropriate 
boundary conditions. The horizontal tiles of the tiling corresponding to the 
diagram in Figure ^ are shown in Figure |21 

It is clear that the positions of horizontal tiles uniquely determine both the 
tiling and the partition tt. The set 

cr(7r) = {(i-j,7ry - (i + j - l)/2)} C Z x iZ (8) 

is precisely the set of the centers of the horizontal tiles. Notice that if (h, t) is 
a center of a tile h + t/2 + l/2\s always an integer. 

Define 

AT AT-l 

i=l i=l 

The image of the inner boundary of our skew plane partitions in the (/i, t)- 
plane is the curve 

h = —B{t ), 

see an example of this curve in Figure^] In particular, the highest layer of hori¬ 
zontal rhombi for an empty plane partitions is the set of points with coordinates 
h = -B{t)-1/2. 
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- 5 - 3-11 3 5 

Figure 3: Horizontal tiles of the tiling corresponding to the partition in Figure 
n]in (t,/i)-coordinates 



Figure 4: Coordinates of corners 
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2.3 Partition function and correlation functions 

Generalizing introduce a probability measure on skew plane partitions by 

Prob({A(t)}) oc (9) 

tez 


where 0 < g* < 1 are parameters. 

We assume that % = 0, t < uq = —a or t > un = b and so the plane 
partition is confined to a a x 6 outer box. The homogeneous case when all 
nonzero qt are equal corresponds to ©• 

For fixed inner shape /i, the partition function is defined by 

2(w.rt= E 

{A(i)}tGZ TT tGZ 


where |7rt| = Y,i^ht+i- 

The correspondence tt < 7 ( 71 ) defined in 0 makes a random skew partition 
a random subset of Z x (Z + i), that is, a random point field on a lattice. This 
motivates the following 

Definition 1. Given a subset U cZx(Z+i), define the corresponding 
correlation function by 

p{U) = Prob(f7 C cr(7r)) = ^ ^ 

7r,t/C(7{7r) ieZ 


These correlation functions depend on parameters qt and on the fixed inner 
shape fj, of skew plane partitions. 

Gonsider the following “local” functions on skew plane partitions : 




1 if {h, t) e cr(7r), 
0 otherwise. 


If U — (^hfi ,)} with tx ^ ^ tji and (^ht ,tx) ^ (.^j fb^ 

correlation function cnj can be written as: 

P{U) = {PhtM ■■■Ph„,t„) = ^'^Pht,tA^) ■ ■ ( 11 ) 

TT teZ 


3 Schur processes 

3.1 General Schur processes 

Schur process, introduced in m is a probability measure on sequences of par¬ 
titions. 
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Parameters of the Schur process are sequences of pairs of functions {(/)f (z)}tgz 
such that is analytic at z = 0 and is analytic at z = oo. For such 

pair of functions <^( 2 )^ consider skew Schur functions 

= det(0+, 

and 

Some basic notions about Schur functions are recalled in the Appendix fXl 
Define the transition weight by the formula 

S'^(A,Ai) = ^Sx/u[(t>^]Sf,/u[(t>~] 

V 

Definition 2. The probabilities of the Schur process are given by 

Prob({A(t)}) = ^ (A(m-i),A(m+i)) , 

mGZ+ 1/2 

where the transition weight is defined above, by (j)t we denoted the pair of 
functions (fy^ and Z is the normalizing factor (partition function) 

Z = n ^)) . 

{A(t)} mGZ+1/2 

If instead of infinite sequences {A(t)} we have finite sequences of length N 
we will say that the Schur process is of length N. 

3.2 Polynomial Schur processes and height distributions 
on skew plane partitions 

We will say that the Schur process is polynomial if functions 4>^{z) are polyno¬ 
mials in . Let us show that the measure 0 is closely related to a polynomial 
Schur process. 

We will parametrize the inner shape as before by assuming that Vi, 1 < f < 
N, Vi G Z are positions of inner corners of the inner shape of the plane partition 
(see Fig. 0 and Ui G (vi, Vi+i), 1 < i < N — 1 are positions of outer corners. 

Theorem 1. The restriction of the measure m to random variables supported 
on subsequences {A(to)}, to < vi, {A(ni)}j 1 < * < A^, {A(to)}, m > vn 

coincides with the polynomial Schur process with parameters 


^m(^) = 

(z - X^) , 

^m(2) = 1. 

(12) 

( 2 ) = 

n 

(2 - , 

(13) 


Vi<m<Ui, mGZ+-^ 


C (Z) = 

n 

(z-x^) , 

(14) 


Ui<.m<.Vi+i , 

mGZ+ ^ 


= 

(z - x +), 

fcW = 1: 

(15) 
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where parameters and qt are related as follows: 


X 


+ 

m+1 


X 


+ 

m 


X^ ^X I 
Ui-^ “i +2 


X 1 a;'*' 2 
— 2 + 2 


a-m+l 



Vi < m < Ui — 1, or m > vn 


(16) 

(17) 


Qvi 1 


(18) 


q 1 ,Ui < m < Vi+i — 1, or m < Vi 


m+ 


(19) 


Proof. Let us restrict the process 0 to the subsequence It is easy to see 

that transition probability from Vi+i to Vi in such subprocess are 

S\\(Vi)^ A(Ui+i) = ^ Sx{yi)/\{ui){^y^J^l|2^ ^iij+3/2’ ■ ■ ■ ’ —1/2) (^0) 

'SA(j;i+i)/A(ui)(a^^._^i_l/2> ^Vi+i-3/2' ' ’ ’ ’ ^«i + l/2) (^^) 


where xff are related to qt as in P|l . Thus this process is a polynomial 
Schur process with (j)f given by II 3 - Conversely, it is clear, that due to the 
identity m any polynomial Schur process can be extended to a probability 
measure on sequences of interlacing partitions with parameters {qt} defined 
as in (tT^ltTHIl . □ 


4 Fermionic representation for correlation func¬ 
tions 

4.1 

For m G Z + i define e(m) = + if Ui < to < Ui and 1 < i < N, and e{m) = — 
if Ui < TO < Vi+i and 0 < i < — 1. This is shown on Fig01 Define = 

{TO|e(TO) = +} and D~ = {TO|e(TO) = —}. 

Let xff be positive numbers related to qt as in (P|l . Notice that for given qt 
the numbers xff are defined up to a transformation x^ x^a^^. 

Theorem 2. 1. The partition function for the height distribution on skew 

plane partitions can be represented as the matrix element of the product 
of vertex operator described in the Amendix VB .S\ as follows 


Z = 


n r_(x+)--- n r+(x-) n r-(a;+) 

UN>Tn'>VN Ui<.m<.Vi^i Vi<.m<.Ui 


n r+(x-)uy\u« 

vi >m>uo 





( 22 ) 
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and 

Z= (1-mi e Z+i 

mi <7712 ,miGD~ ,1712^.0'^ 

2. Assume ti > ■ ■ ■ > then 

{Ph,M ■ ■ ■ Ph„,0 = U n ■ ■ ■ 

7rLGl^-\-^,UM>77l>ti 

ti<m<ti — i ti^i<7n<ti 

n ( 23 ) 

t„>m>«o / 

Here and below ji = hi — B{ti) + 1/2. 

3. Correlation functions m are determinants: 

{Phi.ti ■ ■ ■ Phn.tj/} — detl^K (fi^ hi) {tkt hk)))l<i,k<n (24) 

where 

K{{ti, hi), {t2, / 12 )) = 

1 f f ^-{Z,ti)<^+{w,t2) 

(27ri)2 J\,\<R(t^)JR•it 2 )<\w\ d>+{z,ti)<^_{w,t 2 ) 

^-hi+B(ti)-l/2„^h2-B(t2) + l/2 d'Zdw ^25^ 

z — w zw 

Here |u;| < |z| for h > t 2 , |w| > jz] for C < t 2 , R{t) = min^>t(|a;+1“^) 
and R*{t) = maxm<t(|a;“ |). Functions ^±{z,t) are: 

^+{z,t) = n ( 26 ) 

= n i^-z~^xf,) ( 27 ) 

m<t,7nGD~ ,mGZ+-^ 

Proof. The fact that the partition function and correlation functions for the 
height distribution of plane partitions the matrix element of the product of 
vertex operators as above follows form the formula (Umi for matrix elements of 
products of vertex operators r±(x). 

Using the commutation relations (P)l. and the fact that r_(x)u*^m)o = 0 we 
obtain the product formula for the partition function. 

The operators ifjif* act on the vector as follows: 
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if j = Xi — i + 1/2 for some * = 1,2,... and 

= 0 

otherwise. Using this fact and the formula for the matrix elements of r±(a;) we 
obtain the formula 123) for the correlation functions of densities. 

Moving operators r_ to the right and r+ to the left we obtain the following 
formula for the correlation functions 


{Phi,ti ■ ■ ■ Phn,tn) — 

where 

V'j W = n r_(a;+) Y[ r+(a;/,)"Vj 

n r+(a;-) n r_(x+)-i (28) 

and 

n r_(a:+) n r+(x-)-V* 

n r+(x-) n r_(x+)-i (29) 

Here the operators on the right are given by power series. Commuting formal 
power series gives the following identities: 

r+(a:)“V/cr+(a;) ='ipk - x^k+i 

r_(a:)V’fer_(x)“^ = x^i^k-n 

n>0 

r+(x)-V:r+(x) = y]x-*/L„ 

n>0 

^-{xWkT^ix)-^ =11^1- 

Applying these identities to the formal Fourier transform of tpj (t) and ijj* (t) 
we obtain: 

ij{z,t)= n (1-^2;+)-^ J| (1 - )V>(z) , (30) 

m>t,77iGD'^ m<t,mGD~ 

^*{z,t)= n ( 1 -^ 2 ;+) (1 - z-^a:-)-V*(2;) , (31) 
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where both sides are power series in and are formal Laurent power series in 
z.Uve F, then ^{z)v G z^/‘^F[z-\z]] and 'ip*{z)v = z^/^F[z-\z]]. 

For the inverse Fourier transform of ip{z,t)v and tl){z,t) we obtain the fol¬ 
lowing integral representations: 


“Si/. 


z 'Ip{z)vdz 




(32) 




27Ti , 




^-{w, t) 


^ijj* {w)vdw 


(33) 


Here v G F, both sides are vectors in and these power series converge 

for sufficiently small x^s. 

The contour of integration for z is chosen in such a way that none of the 
poles of ^+{z,t) will be inside of it, this gives \z\ < R{t) = minm>t(|x+|“^). 
The contour of integration for w is such that none of the poles of ^-{w,t) are 
outside the contour. This gives |u;| > R*{t) = maxm<t(|a:“ |). 

As it follows from the Wick’s lemma (|S3 that correlation functions CD 
are determinants of matrices of correlation functions of two Clifford operators 

(EHUSni)- 

{Phi^ti • ■ • Phn,tri^ — 


where 


^a); ^b)) 


(V’ja (ta) '0*, (tb) vo,vo) , a<b, 
~ (^f>) '*^ 0 , '^o) , a> b. 


Notice that a < b iS ta > tb and a > 6 iff to < tf;- Now, substitute (ED and ED 
into Kab and take into account (ED- This proves the formula for correlation 
functions. 

□ 


4.2 

Notice that operators ip{z,t) and 'ip*{z,t) satisfy the difference equations: 

= {I-zx'lj^^^^)ip{z,t), tGD+ 

'^{z, t -p 1) = (1 - z~^x~^^/^)ip{z, t), t G D_ 

0*(z,t- 1) = (1 - za;+_^/2)0*(z,t), t G D+ 

il}*{z, t - 1) = (1 - z~^x~_^j^)-4}{z, t), t G D_ 

These difference equations give the following difference equations for corre¬ 
lation functions: 

Ar((ti, hi), ^ 2 )) — K{(ti — 1, hi -b 1/2), (^2, h2)) 

+ ^q-i/2^((^i “ “ 1/2): (^2,h2)) = F G D+ , (34) 
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hi), {t2, h2)) — K{{ti — \,hi — 1/2), {t2, /i2)) 

+ ~ 1; ^1 + 1/2), (^2, ^ 2 )) = ^ti,t2^hi,h2: ^1 € -D- , (35) 

Using these equations and similar difference equations in t 2 one can express 
all correlation functions in terms of equal time correlation functions. 

Remark 1. The equations m and together with appropriate boundary 
conditions are the equations for the inverse Kasteleyn matrix for the correspond¬ 
ing dimer model. 


4.3 The homogeneous restricted case 

In the homogeneous restricted case 0 < qt = q < 1 ior uq < t < un and qt = 0 
otherwise. 

In the homogeneous case The partition function does not 

depend on a. The functions ^±{z,t) are: 


m>t,m£D+ 

$_(z,t) = n (l-^-V"^a-i) 


In this case we have 

m = 


7n<t,m^D~ 


a ^q Ui < t < Vi 

a~^q~*, Vi-i < t < Ui 


a ^q 
a-^q-^'- 


Ui <t <Vi 
Vi-I <t <Ui 


(36) 

(37) 


R*{t) = 

Notice, that when g —s- 0 the density of horizontal tiles converges to 


p{h,t) 


(27ri)^ 


'|z|=l+e./|u,|=le 


1 dz dw 

Z-W ^h+Bit) + ^^-h-B{t) + i ' 


This integral is 1 when {h, t) is on a “floor” and is 0 when it is on the “wall”. 


5 The thermodynamic limit 

Here we will study the limit g —> 1 of the homogeneous Gauss distribution on 
restricted skew plane partitions when the number of corners in the inner shape 
of diagrams remain finite. 

We assume that q = exp(—r), r —*■ -|-0 and Ui = rui, Vi = rvi and N 
remain fixed and that Uq < Ui < C/i < • • • < Uv < Un- 
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5.1 Asymptotics of the partition function 

It is easy to compute the free energy of the system in this limit: 


F = -\ogZ = - Y. log(l-g”-™) = 


- \\ Y i f log{l - ‘')dfidv\ + o(^) 

Similarly, one can compute the asymptotic of the average volume of a 3D 
partition: 


(kl) = Q 



E 


n — m 

^ _ qn—m 


1 


- i<j ' 




'Vi<u<Ui 


V — ^ 
1 - 



The first formula reflects essentially two dimensional nature of the problem. 
The second formula implies that r~^ is the characteristic length of the system 
when r —> 0. 


5.2 The function S{z) 

Now let us analyze the correlation functions US) in the limit r —> +0. Since r ^ 
is a characteristic scale of the system in this limit we assume t = tir, y = 
remain finite. Depending on the value of (y, r) we will either keep the differences 
ti — tj and hi — hj finite, or we will scale them as appropriate powers of r. 

When r —> +0 the functions in the integral defining correlation functions 
behave as 




= exp 



F{z){l + 0{r)) 


where 


S{z) = [ 
J U 


log(l—z 


)dn- 


fi<r, fiGD- 


} /i>r, 


log(l-ze '^)dfi-{x+B{T))ln{z) 


and F{z) can be computed explicitly. 

In this limit the integral l|23) becomes 


K{hi,ti), (h2, t2)) 


1 

(27ri)^ 




S{z-,XJ-,'ri)-Siw,X2,T2) 

e r 


F{z;xi,ti) y/zw dz dw 
F{w, X 2 , T 2 ) z — w z w 


(1 + 0 ( 1 )). 


(38) 
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The integration contours are described in the previous section. For example, 
if = 2 and Ui < ti,T 2 < V 2 the contours are: 

_ J > \z\,\w\> \z\,\w\> if n > T 2 

z w > \z\ > |w| > e'^ if n > T 2 

The integral can be computed by the steepest descent method. In 
order to do this one should first analyze critical points of S{z) and then deform 
contours of integration accordingly. 


6 Critical points of S{z) and the deformation of 
contours 

The function S{z) can be written as a sum of dilogarithms: 

N N 

S{z) = -(x + B{t)) ln(z) + ^ Li 2 ( 2 ;e“'^‘) - ^ Li 2 (ze“'^‘) - Li 2 {ze~'^) . 
where 


i=0 


Critical points of S{z) are zeros of 


= -{x + B{t)) + 


U 2 {z) = I 








This is equivalent to the following equations. 

/ 7 -)/ \\ 1 1 — zc ,1 — zc , 

-(x+i?(r))-E ^ l^g( 1 _ ) = » 

0<j<i N>j>i 

when Vi < r < Ui and 

/ 7 -)/ w 1 / ^ 1 — 2 : ,1 — zc^^ 

-(x+b(t))- y: i„g(^__j^)-bg(^_^)+ y; iog(__^) = 0 

when Ui < T < Vi+i. 

Exponentiating these equations we obtain 
exp(-x-B(r) -L(t)) ]q {I - ze~^^) = {I - ze~'^) ]q (1 - (39) 

0<i<A^ l<j<N 


where 


_ / s:o<i<*(r7+i Uj) 


if Ei < r < C/g 
+ r - {7j if C/i < r < E^+i 
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The number L{t) is the total length of (—) intervals which are to the left of r. 
It is easy to see that 


L{t) + B{t) = tl2 — uq . 
Therefore the equation can be written as 

eX—/2^ - = f{z) 


where 


/(^) 


aUo 




1 ) 

1 )' 


(40) 


6.1 The number of roots 

Theorem 3. The equation has either N real solutions or N — 2 real solu¬ 
tions and two complex conjugate. 

Proof. The function f{z) has simple poles at z = Vi and z = oo and simple 
zeros at z = Ui. A sample graph of the function f{z) is plotted in Figure |31 



Figure 5: The graph of f{z) for {Uq, Vi, Ui, ¥ 2 , 1 / 2 } = {-2, -0.5,1.5, 2,3} 
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Solutions to the equation are intersection points of the line ^ z — 

gX+T/2 the graph of the function f{z). Any straight line of non-infinite 
slope obviously intersects the graph of f{z) in at least N — 2 points (and in at 
most N points, since the degree of / equals N). Hence among the roots of 14011 
there is at most one complex conjugate pair. □ 

6.2 The N = 1 case 

First, consider the equation When N = 1 the normalization of U’s and H’s 
given by the equation implies Vi = 0. Introduce variables 

X = T = e"^, a = , /? = 

The equation for critical points of S{z) is quadratic: 

iP - XT)z^ + {X + XT-{1 + aP))z -h (a - AT) = 0 (41) 

The discriminant of this equation is: 

X = {X-XT-a + Pf -{X + XT)2{1 - a)(l - /3) -b (1 - a^)(l - 0^). (42) 

The structure of solutions depend on the value of the discriminant A. 

• When A < 0 there are two complex conjugated critical points 

• When A > 0 these critical points are real. 

• When A = 0 the two simple critical points degenerate into one double 
critical point. 

6.3 The N = 2 case 

This is the smallest value of N when the function S{z) can have a triple critical 
point. 

The function S{z) for the case when N = 2, and Hi -b V 2 < t < V 2 is 

pVi pr nU 2 

S{z)= \n{l-z-00dfi+ z-00dfj,- ln{l-z-00dn-{x+B{T))\n z 

JUq Jvi-\-V2 Jv2 

We choose branches of logarithms such that the derivative of S{z) has brunch 
cuts e"^], and [e^^,e^^]. 

The function S{z) has three critical points. They are either all real or there is 
a complex conjugate pair of simple complex critical points. Geometrically these 
critical points correspond to intersection points of the line —with 

the graph of the function f{z). 

When the line z — intersects the graph of the function f{z) 

transversally, the intersection points are simple critical points of S{z). 

At double critical points the the line z — is tangent to f{z) but 

it does not bisect the graph of f{z) at the point where it is tangent to the graph. 
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At a triple critical point the line is tangent to the graph of f{z) and bisects 
it. 

By the definition, a double critical point z of the function S{z) satisfies the 
two equations S'{z) = 0 and S''{z) = 0. We have 

d 

z —^(z) = -x-ln(ze"'5 -e^) + ln/(z) 




ze~i — f{z) 


This gives the system of equations for double critical points: 


eX-rt^z- ^X+r/2 ^ j^z) 
gXo-ro/2 ^ j,^z) 


(43) 

(44) 


These equations define the curve in the (x, r) plane. We will say the point 
(x, t) on this curve is generic if it is not a triple critical point, i. e. if S"'{z) ^ 0 
where z is the corresponding double critical point of S{z) (a solution to (HHll . 

Denote by zq the triple critical point of S{z) and by (xo; T)) the correspond¬ 
ing values of x and r. By definition 5''(zo) = S''{zo) = S'''{zo) = 0, which gives 
one more equation in addition to 63 : 


f'izo) = 0 


It is clear from the shape graph of f{z) that for each value of r there are either 
2 or 4 double critical points of S{z). They satisfy the following inequalities: 

1. —oo < T <Vi, then zi < 0 and 0 < Z 2 < e'^^ 

2 . Vi < T <Ui, then zi < 0 and < Z 2 < 

3. C/i < r < To, then zi < 0, e"^^ < zi < zq, zq < Z 3 < and 
e"^ < Z4 < 

4. To < r < V 2 , then zi < 0 and < Z 2 < eX'^ ■ 

5. Vi < T <Ui, then zi < 0 and eV'^ < Z 2 . 

It is also clear that if tq > C/i, then 


6.4 Deformation of integration contours 

We want to deform contours of integration in 63 to position them in the way 
the steepest descent methods requires. If Zc is a critical point of S(z\ and if it 
does not lie on a branch cut of the function 5'(z), the integration contour should 
be deformed to a contour which lies on the curve 5(5'(z)) = Q(S'(Zc)). 
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The function S{z) has branch cuts along the real line. However, only 
exp(^^^)F(z), which is the leading term of the asymptotic of 
appear in the integral (P|l. 

Zeros of $+( 2 ) are accumulating along Therefore we can not 

deform Cz through this segment but we can deform through any other part of 
the real line. 

Similarly, zeros of <&_(?«) are accumulating in the segments and 

(e’^be'^). Thus, the contour can not be deformed through these segments 
but can be deformed through any other segment of the real line. 

Therefore we have to deform contours Cz and to the union of appropriate 
branches of curves 9(S'(z)) = 3(5'(zc+*0)) and9(S'(z)) = 3(S'(2:c—*0)). Figures 
nci and IHl show these curves for simple critical points, double critical points 
and the triple critical point, respectively. 



Figure 6: Level curves of 3(5) when all critical points are simple 



Figure 7: Level curves of 3(5) when there is a double critical point 
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Figure 8: Level curves of 3(5') when there is a triple critical point 


Deformed contours of integration and are shown in F igures El IHl 

nnmnd for simple, double, and triple critical points respectively. 


Figure 9: Deformation of the integration contour Cz for a double critical point 

While deforming contours Cz and one should keep track on the residues 
at the pole z = w. We will discuss this later. 

7 Simple critical points and bulk limit 

7.1 Bulk limit of correlation functions 

Here we will compute the asymptotic of the integral when there is a pair 
of complex conjugate critical points and when Ah = hi — h 2 and At = ti — t 2 
are fixed in the limit r —> 0. 
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Figure 14: Deformation of the integration contours for a simple critical point 
when ti <t 2 


Deform contours C'z and to Ci and C 2 respectively as it is shown on Fig. 
m for ti > t 2 and as it is shown on Fig. m for ti <t 2 . 

Taking into account the residue at z = w we have the identity: 



Here we take + for At > 0 and — otherwise, li Ui < t < Vi+i then 7 + is a 
simple curve connecting Zc and Zc and passing through positive part of the real 
line with < z. Similarly, 7 _ is a simple curve connecting Zc and Zc and 

passing through the negative part of the real line, li Vi < t < Ui then 7 ± are 
again simple curves connecting Zc and Zc in such a way that 7 _ intersects the 
negative part of the real line and 7 + intersects the positive part of the real line 
at < z. 

As r ^ 0 only the second term in the right hand side of 103) will have finite 
limit. The first term will vanish. The computations are similar to |15|. The 
integrals along 7 ± converge as r —> 0 to 


K{Ah, At) 


Biy^}"'^\At, Ah + e{T)At/2) , At >0 

Bj^}'^'^\At, Ah + s{T)At/2) , At < 0 


(46) 


where 

' 27rz 

Here £(t) = 1 when Vi < t < Ui and £(t) = — 1 when Ui < t < Vi+i. The 
contours 7 ± are as above. 

For the limit of correlation functions we obtain the following answer: 


lim {phi tl • ■ ■ Phn,tn) 

r^O ~ ~ 
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where 

^a,h — ^((ha {td — lim K{^{tdi ^a)) 

r—»-+0 

Correlation functions for finite q satisfy recurrense equations and C^5II . 
In the limit these recurrence relations turn into difference equations for corre¬ 
lation functions CD: 

K(h, t) - K{h + £(t)1/2, t - 1) -I- - £{r)l/2, t - 1) = Sh^Stfi 


where e(T) = -|-1 when r G £)+, e(T) = — 1 when t £ D . These equations can 
also be directly deduced from the integral representation of correlation functions. 

Remark 2. The difference operator in these equations is the Kasteleyn matrix 
on the infinite hexagonal lattice. The pairwise correlation function EH) can be 
regarded as the inverse for this matrix with the boundary conditions determined 
by (t,x). 

For one-time correlation functions we have: 


K{h, 0 ) 


I —h 


sin{9h) 

nh 


where 6 is the argument of Zc, that is, Zc = \zc\e^^■ 

Notice that the factor \zc\^ does not contribute to the equal-time density 
correlation functions and we have: 


lim ■ ■■Ph„,t) = det 

r— 

Here we assume that hi > /12 > ... /i„. 


/ sin(6>(ha - hb)) \ 

V n{ha-hb) )i<a,b<n 


7.2 Limit shape 

The form of the one-point correlation function implies that the density of hori¬ 
zontal tiles is 

P{r, x) = - 

TT 

where 9 is the argument of Zc- The limit shape can reconstructed from this 
density by integration: 

z{t,x)= [ (l-P(Ts))ds. (47) 

J-b(T) 

x{t, x) = z{r, x) - X - ^ , y{r, x) = z{r, x) - X + ^ , (48) 

Here &(t) = max(— t/2 -|- Uq, t 12 -|- U-^). Thus, the information about the limit 
shape is in the structure of critical points of the function S{z). 

When the point (r, y) is such that all critical points of >S'(z) are real the limit 
of correlation functions EH is either 1 or 0. It is zero if the maximum of S{z) 
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is inside of the cycle of integration and it is 1 if it is outside. The corresponding 
point {x,y,z) lies on a facet (flat part of the limit shape). 

If (r, x) is such that there is a pair of complex conjugate simple critical 
points, the point {x, y, z) lies in the disordered region (curved part of the limit 
shape). 

The frozen boundary (that is, the boundary between the disordered region 
and the facets) limit shape correspond to (%, r) for which there exists a real zero 
of S{z) of the multiplicity at least 2. Cusps of the frozen boundary correspond 
to triple critical points. 

We plotted the frozen boundary in the (t, x)-plane on Fig. ^jfor N = \ 
and on Fig. ^|for N = 1. The vicinity of the cusp is magnified on Fig. El It 
is instructive to compare these curves with the result of numeric simulation in 
Figure 121 



Figure 15: An example of frozen boundary for N = 1 



Figure 16: An example of frozen boundary for N = 2 
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Figure 17: The vicinity of a cusp 


8 Double critical points and the scaling limit 
near the boundary 


8.1 


Now, let us assume that (xoih)) are such that zq is a double critical point of 
S(z). Consider the vicinity of (zo,Xo,to) with coordinates z = zoexp(^),x = 
Xo + = ro + St. Lowest degree terms in the Taylor expansion around 

(zo,Xo,To) are: 


S(z,x,r) = S(zo,xo,ro) + - (Sx-BS t)^ + ^CSt^^ + ^DSt^^ 

+ —ESt^ + HSt + -^GSt"^ — ln(zo)dx + ... 


where 


oz f{zo) 

a2 _ ZQ. , ZQ. 

n _ ^ M _ 2oe = + e 2 

B — z „ „ d(z)|2q^to — 


drdz 


2(zoe ° — 62 °) 


zo 


^ (zoe-^-e^)2 


Zq 


d d'^ 

^ = ~zoe-^ - e^ f 


qZ 

E=^,S(^z% 


Zo 


(zoe 2° —6 2° )2 


(49) 

(50) 

(51) 

(52) 

(53) 


E — Q^S{z)\zo,To 


-i In(zo) + ln(zo - tq £ 

-iln(zo) + ln(e^ - Zo) - To, t-q G T>+ 


(54) 
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.•^0 


To e £>_ 
To e D+ 


(55) 


G* — '^(^) I 20,To 


zo-e^o ’ 

zq 

€■^0 —Zq ’ 


Notice that C = —D, E = —D, D > 0, A/B > 0. The sign of A depends on the 
nature of the interface. It is positive if the frozen region is above the melted 
region and it is negative otherwise. It changes signs at the points zq = ex.p{Ui) 
and at the triple critical points where f''{zo) = 0. 

Rescaling local coordinates 6x and 6t as 


^ = r^a, Sx — BST = r3x, ST = r3y 


we have 


S{z,x,t) - S{zo,xo,tq) 
r 


-xa+ ^Cya'^ + ^Dy'^cr + ^Ey^ + o(l) 


8.2 

As r ^ 0 the leading asymptotic of the integral is determined by the leading 
asymptotic of the integral: 


K{hi,ti), (/i2,l2)) = 


(27r*)^ 


S(2;xi ,zi )-S(u,,X3,z,) ^/zw dz dw . I/O, 

e -(l+0(r^/'^) 


'C^ JCn, 


z — w z w 


If Vi < T < V 2 , the contour of integration is 


(56) 


f 1 > |z| > Iwl > ii Vi < T < Ui 

\ > \z\ > |w| >6"^ if t/i < T < V 2 


for ti > ^ 2 - When ti < t 2 the only difference is that \z\ < |u>|. 

Assume that coordinates {hi,ti) are in the vicinity of the point (xojTo), and 
that as r —> 0 they scale in the following way: 

- i i 

Xi = 'rhi = Xo + r^Xi - Br^yi, n = rU =ro-\-r^y^. (57) 


In this limit condition ti >^ 2,^2 > ti translate to yi > j /2 and j /2 > Vi 
respectively. 

Deform contours of integration in such a way that they will pass through 
this critical point and will follow the branches of curves ^(^'(z)) = 5(5'(zo)) 
and ^{S{w)) = Q{S{zo)) where ^^(^(z)) and —^{S{w)) have local maximum. 
Deformed contours are described in the section EH 

The leading contribution to the asymptotic comes the vicinity of the point 
zo where we will use local coordinates z = zo{l + ra a) and w = zo{l + ran). 
The saddle point integration contours for a and k in the limit r ^ 0 are shown 
on Fig. IT^and Fig. ^jfor A > 0. 
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The leading term of the asymptotic of the integral iCT) is given by 


K{{hi,ti), (h2M)) = exp(^(2/i - 2/2) + ^(l/i - 2/2) - - 3^2) 

-—^^iyi-y^) + J^^Eiyi-V2))KA{{.xi-^yl, yi), {x2-^yl, y2)){l+0{r^/^)) 

( 58 ) 

For yi > 7/2 the function Ka is 


KA{{ai,hi), ( 02 , 52 )) — 
1 


f f A D 

/ exp( —(cr^ - - —(6icr^ - 62^^) - aio-+ 02^) 

0 ic+ ic- 3 ! 2 


dadn 


a — K 


( 27 ri)- jc<+ 

For yi < 7/2 it has an extra term which comes from the residue at cr = r: 


( 59 ) 


1 f f A D 

KA{{ai,bi),{a2,b2)) = / exp(—(cr^ - k^) - —(fepT^ - 62^^)- 

(2717)^ Jc+ Jc- 3! 2 

, dadn 1 /" r n 1 \ 2 / \ \ 1 rnr^\ 

oiCT + a2/c)- 'r TT- / exp(—— (Oi — 02)cr — (oi — a2)cr)acr. ( 60 ) 

a — K Ztti 7c+ ^ 

The function i 4 r^((xi, 7/1), {x2,y2)) is discontinuous at 7/1 = 7/2- The discon¬ 
tinuity is determined by the second term. It is not difficult to see that it is equal 
to 


1 


27 rD|yi - y 2 


{Xi - X2f . 
2-D|yi - y2| 


Notice that since correlation functions 123 ) are given by determinants, the 
exponential factors in front of the integral in 1(23 will be canceled and we have: 

lim {ph^M ■ ■ ■ Ph^,u) = dei{KA{{Xa - ^yl, ya), {xb - yt, yb))) 

r —>0 2, I 

Here we assume that {hi,ti) are scaled as in |23 
It is easy to verify that 


KAi{xi, 0 ), {X 2 , 0 )) = 


Ai(Aa:i) Ai^(Aa:2) — Ai(Aa;2) Ai^(Axi) 


where Ai(a;) is the Airy function: 


Ai(a:) = 




A(a;i - X2) 


exp(-t^ -|- ixt)dt, 

> «J 


and A = 
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For the scaling limit of the density of horizontal tiles we have: 


PA{x,y) = 
1 


— i 


(27t)2X 


exp —i 


R JR 




3! 


AV3 V 2A 


— a; I (s — f) 


dsdt 


s — t — iO 


( 61 ) 


It can also be written as 


PA{.x,y) = 




Ai{v)^dv 


In this form it is clear that the density function is positive. 


9 Triple critical points and the scaling limit near 
the cusp 

9.1 The singularity of the limit curve 

Then the following holds if Zq is a triple critical point 


_5 

'dz^ 


{z^)^S{z)U^=z^o 


f{zo) 


(62) 


The proof is straightforward. First, one computes the fourth derivative: 


{z§-/S{z) = — 
oz [ze 


2z^ 


r 






2 — 62 )2 


{ze 2 — 62 )3 


— + {z2^ + 3z^^y + zH^r) 


/ 


/ 


/ 


This expression reduces to (inu after taking into account equations for zq- The 
second derivative of f{z) vanishes at zq and, as it follows from the graph of f{z), 
it is negative for z < zq and positive for z > zq. This implies the positivity of 
/(^^(zo). This means that the sign of A is determined by the sign of f{zo). The 
later is negative if zq < exp({7i) ( when the cusp on the limit shape is turned 
to the right) and is positive otherwise (i.e. when the cusp is turned to the left). 

Now, let us find the behavior of the boundary curve near the cusp. Assume 
that Zq is a triple critical point corresponding to the cusp with singularity at 
(XOiA))) be. zo,ao and bo satisfy the system: 


gXo-ro /2 ^ 

(63) 

gXo+To /2 ^ Zof'{zo) - f{zo), 

(64) 

II 

o 

(65) 


If a point (x, t) is at the boundary curve, it satisfies the equations: 


gX-r/2 ^ gX+r/2 ^ 


( 66 ) 
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Let (z, x, t) be a double critical point in a small vicinity of (zq, xo, to): 


X = Xo + Sx, t = to + St, z = zo+e 

Then from the equations we obtain the following asymptotic of the 

boundary curve near the cusp: 

exp(xo - To/2){6x - 5 t/2) = + 0{£^) 

2 D 

exp(xo + To/2){Sx + St/ 2) = ^Zof^^'> (zo)£^ + ^ (zq) + 2f^^'> {zo))£^ + 0{£‘^) 

2 D 

when 6^0. From here we have: 

Sx = ^exp(-xo - y)((zo + exp(ro))/(^H^o)e^+ 

i((zo + exp(To))/^^H^o) + 2/(3)(zo))e^ + ■•■) (67) 

i exp(-xo - y)((^o - exp(ro))/(^)(zo)e^+ 

i((zo - exp(To))/^^Hzo) - 2f^^\zo))£^ + ■■■) (68) 

After reparametrization 

(zo-e-°)/W(zo)-2/(3)(^o) 

6/(3)(zo)(2:o - e^“) 

we have: 

St = iexp(-xo - y)(^o - exp(ro))/^^Hzo)CT^ + O(cr^) , 

Sx= ^exp(-xo-y)((zo+exp(To))/(^)(zo)o-^-^^^^^^h5^/(3)(^p)cr^)+0(cr4) . 

This is a parametrization of a cusp singularity in the boundary of the limit 
shape. 

9.2 The asymptotic of correlation functions 

Expanding the function S(z, x, t) near the triple critical point zq we obtain the 
following lowest degree terms of the Taylor expansion: 

S{z, Xo + Sx, To + St) = S{zo, Xo, ^o) + 2^'^ _ Sx^ + BSt^ + ]^CSt^'^ + 

+ 2 ^St^ — ln(zo)^X + BSt + o(l) 
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where 


A = 


f{zo) 


and B and C are as before, given by Recall that A < 0 is the cusp in 

the limit shape is turned right and ^4 > 0 if it is turned left. 

Scaling local coordinates as 




Sx — BST = r^x, 6T = r'^y 


we obtain the following asymptotic for S: 


S{z,x,t) - S{zo,xo,To) 
r 


A 4 C 2 D 2 

= -cr -xa+jya + jV + 
H ^ In(zo)^ , 

^1/2 j,l/4 ^ 


B\n{zo) 

7*3/4 


y + o{l) 


Now let us find the asymptotic of the integral (ISSt as r —!■ 0 assuming that 
coordinates are scaled as 


3 i i 

= Xo + - Br^yi, rU = tq + . (69) 

where (xoj'b)) are coordinates of the tip of the cusp. 

The asymptotic of the integral itCT) as r ^ 0 and 

rhi = Xo + r^x^ - Br^yt, rU = tq + r^yi ■ 

can be evaluated by the steepest descent method. It is determined by the 
contribution from the triple critical point zg and after deforming contours of 
integration as it described in section the leading term of the asymptotic is 
given by the integral 


K{{hi,ti){h2,t2)) = exp{-^-^^^{xi-x2)+^iyi-y2)+ ^^^J}l°\ yi-y2)+ 

^{yl - yl))Kp{{xi,yi), {X2,y2)){^ + (70) 

where for yi > y 2 we have 
Kp{{xi,yi),{x2,y2)) = 

1 f f 4 nC', 2 2 n .dadn 

.O ■N 2 / / exp — cr -K) + —{yi(T - y2K ) - Xia + X2K) -h 

[z-my Jc^ Jc2 2 cr - K 

\j 2 TTD\yi - y2\ 2 D\yi - y2\ 

For yi < j /2 there is an extra term coming from the residue at cr = r: 
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Kp{{xi,yi), {x2,y2)) = 

1 f f C, ^ o, .dadn .. 

fr) ■N2 / / exp — cr^-K^ +—(2/icr^-?/2K -a;i(T + a;2K)-. ( 72 ) 

{2Tny Jc^ Jc 2 4! 2 cr - k 

Integration contours for > 0 are shown on Fig. 1211 and 1201 


/\ 


Figure 20: Integration contour for k 



Notice that the exponential factors cancels in the limit of correlation func¬ 
tions and we have: 

lim, {Phuti ■ ■ ■ Ph„,tJ = det{Kp{{xa, ya), (Xb, Vb))) 

r —*^0 

Here we assume {hi,ti) scale as inlHIHl. 
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10 Some properties of Pearcey kernel 

Recall that the asymptotic of correlation functions is determined by the follow¬ 
ing integral, which we will call Pearcey kernel 


Kp{{xi,yi),{x2,y2)) = 

1 f f 4nC'/9 9n .dadn . 

.0 02 / / exp(— CT -K) + —{yia - y2K ) - xia + X2 k) -. ( 73 ) 

{2Tny Jc^ Jc2 2 a-K 

After the appropriate scaling of variables we can set A = 6 and (7=1. 
Below we will review some of its properties. 


10.1 

The function Kp{{xi,yi), {x 2 ,y 2 ))) (with A = 6 and C = l)satisfies the follow¬ 
ing differential equations: 


for f = 1,2 and 


where 


d 1 


(^ + = -P+{xi^yi)P-{x2,y2) 


1 f 1 

P±ix,y) = ^ exp(±(— -f -j/cr^ - xa))do 

For functions P± we have: 


d 

i^±^ + y-)P^=0 
d 1 


10.2 

The function Pp(x,y) can be written as 


P+{x,y) 





) 


or, as 


P+{x,y) 



)( 


2\( ^xuis 


- e-^^yds) 
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From this integral representation one can find a power series formula for the 
integral: 


P+ix,y) = 


E 


i>o,m>o,n+m=0(2) 


/ 1 \ on-l-p/n+m+1 ^ 

(2n + ly.ml 


or, 


p+ix y)= ^ V (-l)^+^2"^r(fc + I + 1/2) ^ 

+ ^ (4A: + 1)!(20! ^ ^ 

j_ (_i)fc+^2^^+^r(fc + l + 3/2) 4fc_)_3 2;+i 

’ ,h (« + 3)!(2i + l)! " 


Now, using the identity r(A: + l/2) = 2 arrive to the formula 

-21 

{Ak + l)\{2l)\{k+1)\ 


P,(xy)= ^ (-l)fc+^2 ^^+i(2fc + 20! 


k,l>0 


2 


Y ^ (- l )*^+^2 ^^( 2 fc + 2 ? + 2 )! 4 fc _|_3 2 ;+i 

/ > //I 7, I 0M/O7 I lM/7, I 7 I IM*^ ^ 


k,l>0 


(4A: + 3)!(2Z + l)!(fc + Z + l)!‘ 


10.3 

The integral 


1 1 

P-{x,y) = ^- J + -jyt^ + ixt)dt 


can be easily expanded into a power series in x and y. 

„ , ^ 1 V- (-l)"2"-3/2r(H^ + 1/4) „ 

P-{x,y) = - ^ -X y 


n,m>0 


(2n)!TO! 


10.4 

One-time correlation functions can be expressed in terms of P± as follows:: 

K{xi,X2\y) = 

Pl{xi,y)P-{x 2 :V) - P+{xi,y)P'_{x 2 ,y) + P+{xi,y)P”{x 2 ,y) + yP+{xi,y)P-{x 2 ,y) 


Xi — X2 


(74) 


Here P'{x,y) = -^P{x,y). This identity follows from 

exp(Ecr'‘ - k"^) + Eyicr^ - y2H^) - xia + X2 k) = 0. (75) 
4! 2 a — K 
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Taking the limit xi ^ X 2 = x m. iTTl) we obtain the following expression for 
the scaling limit of the density of tiles near the cusp: 

p{x,y) = P\_{x,y)P”{x,y) - P'l{x,y)P'_{x,y) - xP+{x,y)P_{x,y) 

This function is plotted on FigEa 



Figure 22: The scaling limit of the density of horizontal tiles near the cups 


A Schur functions 

Recall that a partition is a sequence of integers 

A = (Ai > A 2 > • ■ • > Aji > 0) 

A diagram of a partition A (also known as its Young diagram) has Ai boxes in 
the upper row, A 2 boxes in the next row etc., see Figure na Unless this leads 
to a confusion, we identify partitions with their diagrams. The sum 

|A| = Ea 

2>1 
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is the area of the diagram A. 


Figure 23: A Young diagram 


A skew diagram \/^ is a pair of two partitions A and /i such that Ai > 
Ml) ^2 > M 2 ) • ■ • ■ Graphically, it is obtained from A by removing mi first boxes 
from the first row, M 2 first boxes from the second row etc. The size of the skew 
diagram A/m is 

|A/m| = |A|-|m| 

A semistandard tableau of the shape A/m with entries 1,2,3,... is the result 
of writing numbers 1, 2,3,... in boxes of the diagram, one in each box, in such 
a way that the numbers weakly descrease along the rows and strictly descrease 
along the colums, see an example in Figure !^ 


1 1 1 1 2223] 

2 2 3 3 3 [ 
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5 5 6 



Figure 24: A semi-standard tableau 


The Schur function corresponding to the skew tableau A/m is a symmetric 
function of variables xi, a; 2 ) ■ ■ • which can be defined as the sum over all semi¬ 
standard tableaux of shape A of monomials in xm 

SX/t,{xi,X2, ■■■) = ^ X^^X™"" . . . 

'Xx/ fj. 

where rrii is the number of occurrences of i in the tableau. 

A semi-standard tableaux of shape A/m can be identihed with sequences of 
skew diagrams in the following way. 

Let us say that A m if M interlace, that is, 

Ai > Ml — -^2 ^ M2 ^ ^ • 

Let us call the sequence of diagrams A(l), A(2),... A(n) interlacing if A(l) A 2 
■ ■ ■ >■ \{n). It is clear that if we will associate with such sequence the semi¬ 
standard tableau with |A(1)| — |A(2)| entrees of n, |A(2)| — |A(3)| entries of n— 1, 
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etc. It is also clear that this correspondence gives a bijection between semi- 
standard tableaux of shape A(l)/A(n) and interlacing sequences of diagrams 
which start with A(l) and ends with A(n). 

For skew Schur functions this bijection gives the following formula 


• ■ ■ ; 0 , 0 , ... ) — 

E 


JA|-|A(1)LIM1)|-|M2)I |A(n-l)|-|/x| 


(76) 


A>“ A(l) >“ •••)>-A (n —1) ^ 


from now on we will write SA/;i(a:i.a; 2 ,...,x„) for sa(xi, ..., a;„, 0, 0,...). 
Notice that 


SA/^(xi,o,...) = ri 


^ ^ A, 


(77) 


and therefore we can write 


^A//i(^l5 ■ ■ ■ 5 ^n) — 

E ■SA/A(l)(a^l)sA(l)/A(2)(a^2) • ■ ■ SA(n-l)//i(2^n) (78) 

A;^ A(l)>->-A(n— 

The function sx{xi,X 2 ,... ,Xn) is the character of the irreducible representa¬ 
tion of GLn with the highest weight A computed on the diagonal element with 
entries xi,X 2 , ■ ■ ■ ,Xn- The formula IIH) is the result of the computation this 
character in the Gelfand-Tsetlin basis. 


B Semiinfinite forms and vertex operators 

B.l Semiinfinite forms 

Let the space V be spanned by fc, fc G Z -|- i. The space F = V is, by 
definition, spanned by vectors 


US = SiAS2AS3A... , 

where S = {si > S 2 > ■ ■ ■} C Z + ^ is such a subset that both sets 

5+ = 5\(Z<o-i) , ^_ = (Z<o-i)\5 

are finite. We equip A“ 2 “ V with the inner product (.,.) in which the basis {us} 
is orthonormal. The space F is also called the fermionic Fock space. 

The infinite Clifford algebra Cl{V) is generated by elements k G 

Z -|- i with defining relations 

ipkipi + = 0 , tpltpi + tplipl = 0 , = 5k,h 

It acts on the Fock space F as 
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4’k ^£1A £2 A ... S£ A fc A Si+i A . 


■) = 


(-l)'si A £2 A . ..^Asi+i A ..., 


tpl (£]_ A£2A£3A...) =0, fcG Z\S 

The space F is an irreducible representation of Cl{V). Notice that the 
operator representing ipl. is conjugate to the operator representing ipk with the 
respect to the scalar product in F. 

The Lie algebra of Z x Z-matrices with finitely many entries acts natu¬ 
rally on V and therefore acts diagonally on semi-infinite wedge space F. This 
action extends to the action of Ooc El and is reducible. Irreducible components 
F{m) are eigenspaces of the operator 


k<0 k<0 

The Fock space decomposes into the direct sum 


F = ®m&F{m) 

of irreducible representations of a^o m 

The subspace F(to) is spanned by vectors S 1 AS 2 AS 3 A... with =to— i-|-l/2 
for sufficiently large i. It is generated by the action of Ooo on the vacuum vector 
( the ttoo highest weight vector in ): 

= to — 1/2 A TO — 3/2 A TO — 5/2 A ... 

The operators ipk and ipl act on as 

= 0, A: < TO - 1/2, (79) 

= 0, A: > TO-1/2, (80) 

They can be regarded as Ooo-intertwining operators if) : V ® F^™^ ^ 

where V = C ^~^2 is the vector representation of 

Vectors in the space F^"*^ can be parameterized by partitions. For a partition 
A define 

=Ai-|-to— I/ 2 AA 2 +TO — 3/2A... 

It is clear that these vectors span F^™). 


39 









B.2 Clifford algebra and vertex operators 


Consider elements 


^ ^ '4^k 

They satisfy commutation relations 

[o^rij Otrn\ — 

[an,'0fc] = 

[a„,V'fc] =- 


n = ±I,±2,... , 


77- Sn,—mj 

(81) 


(82) 

i^k-n 

(83) 


It is clear that 

= 0 

for n < 0 and m G Z. Vertex operators are the formal power series 


r+(x) = exp(^ r_(x) = exp(^ 


a>l 


a>l 


X 

n 


The operator r_ (x) acts finitely in the space F ( applied to any vector of F it 
acts as a polynomial in x). In particular 

r-WV’=V‘’ 

The operator r+(a:) is conjugate to r_(a:): 

(r_(a:)f,w) = (f, (84) 


and since the scalar product is symmetric 

(r_(a:)t',w) = (r+(a:)w,w) 

Notice that its action is dehned in F not only as a formal power series in x. 
In a weak sense operators F ± are operator- valued functions which are analytic 
at a: = 0. 

Define the formal Fourier transform of i/'fe, V’fc Eis power series 


V’(z) = Yl V'*(2) = Y ^ (^^) 

fcez-i-i/2 fcez-i-i/2 

These operators and vertex operators satisfy the following commutation re¬ 
lations: 

r+(a:)r_(2/) = (I - a:2/)r_(2/)r+(a:), (86) 

T+{x)'ip{z) = (I- z~'^x)~'^tp{z)T+{x) (87) 

r-{x)'ijj{z) = {I — xz)~^'tp{z)T-{x) (88) 

T+{x)'tp*{z) = {1- z-^x)i;*{z)r+{x) (89) 

r-{x)'>p*{z) = {1 — xz)'ip* (z)r-{x) (90) 

Here left and right sides are corresponding formal power series. 

The following is a well known statement. 
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Theorem 4. The following identity holds: 

(r_(a:i)r_(a;2)... r_ = sx/f,{xi,.. .,x„) ( 91 ) 

This is a well known statement, for a proof see for example HU. we will give 
a proof of it here. The key step is to show that the identity EU holds for n = 1 
which can be easily derived from the fact that 

where ji = Ai + to — 1/2, j2 = A2 + m — 3/2,.... and from the identity 


r_(x)ip(zi).. = ]^(1 - a;Zi) ^ip(zi)...if(z„)v^’^ 

i=l 


for generating functions. 

For matrix elements of generating functions 'ip{z) and if*{z) we have: 




C’) = ar*''" 
-r’) =(s)’"'"' 


1 

1 — zlw 
1 

1—w/z 


, |z| < |w| 

, |z| > |w|. 


(92) 


These identities follows from the summation of a geometric series. 

Let and Bf = YljBfj'if* for a = l,...,n. The following 

identity is known as a Wick’s lemma: 


where 


[AlBlAlBl... AlBlvt\vt^) = det(K,,) 


l<a.6<ri 


/c, = 




a < b 
a > b 


(93) 


C Some asymptotic for limit shapes 
C.l 

The frozen boundary is singular at t = Vj. When t —>■ Vj the singular branch 
of the boundary curve behave as 


It - V, 


i-i 


X(r) = 2ln(^^) + iy i/2 + ^(y+i - U.)+ 


2=0 


, ,Y{kA3-l<k<N\^ 6 + 

tT ~ )^0{\r-V,\) . 


n 


l<k<N 
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C.2 


The limit curve is tangent to the lines t = Uq and r = Upf at the points (xl, Uq) 
and {xr, Un) where 


XL 


N 


i:i"( 

i=i 


1 _ Uo 

l_e^o-u,)+ 2 


^ 1 _ 

i=i 


2 


Ifr = ?7o + eorr = 17Ar — e and e ^ +0 the asymptotic of two double critical 
points is 


r e'^o(l ± + 0(e)) if T = f7o + e 

\ e^^(l ± + 0(e)) if T = C7^ - e 


where Cl and Cr are some constants. 

The boundary curve near this point behave as: 

= / Xl(1 + \ /{t - UqD l + 0{\Uo - r)|, Ht^Uq + O 
1 Xi(l + V {Un - tDr + 0{\Un - t|), if t ^ t/AT - 0 

where, again, Dr and Dr are some constants. 

Similar solution exists near each point T = Uj,j = l^...,N—1. 


C.3 

Let us find the asymptotic of the density function near the left boundary of the 
limit shape. Assume 

X = Xl + 5x, t = Uo + 5t 

for some positive Jr ^ 0. 

Solutions to iM) have the asymptotic z = exp(C/o)(l + 5z). Let us find 5z 
as a function of Jy Jr. We have the following asymptotical expansions: 


N 


n 


1 - ze 
1 — ze~^i 


N 


n 


1 - 

\ — ^Uq — Uj 


(1 - C5z + 0(Jz^)) 


where 


N 




eUo-Vj 

X — ^Q — Uj 


eUo-Uj 

X — ^Uo — Uj 


Keeping leading orders in <5%, Sr, and Sz in (EIHi we obtain the equation for Jz: 


CSz^ + Jz((C + l/2)Jr — Jx) + Jr = 0 
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In the region 5t cx (^x)^ have two asymptotical solutions 


5x 


From here we obtain the asymptotic of the density function in this region: 


ISx^ 

St 

/ 2 

c 


p = — ~ arctan(4 


' 4:6t 


- 1 ) 


Here 9 is the argument of 5zi. Notice that p 1/2 — 0 as dy ^ +0- 


D The symmetry of correlation functions 

Change variables in the integral 
hi), (t2, /i2)) = 

1 f f d’_(z,fi)d>+(w,f 2 ) yiw dzdw 

(27ri)2 jR^(t„)<\w\<l ^+{z,ti)^_(w,t 2 ) z-w zw 

(94) 


from z to w ^ and from w to z It becomes 
K{{ti,hi),{t2,h2)) = 

1 f f $-(z, -t 2 )^+iw,ti) yiw dzdw 

J\w\-^<rmn{l,R{tl) (-2, ^ 2 )$-(w, H ) Z-W ZW 

(95) 


where 

^+iz,t) = ]^(1 - 2i+) 

m 

z~x;;,) 

m 

and = xtm- 

Thus, we have the following “reflection” symmetry of correlation functions: 

^((jl,H),(j2,t2)) = K{{j2,-ti){ji,-ti)) , 

This symmetry is obvious on the “microscopical level”. It corresponds to 
the reflection of the tiling in t-direction. 


43 



References 


[1] M. Aganagic, A. Klemm, M. Marino, and C. Vafa, The Topological Vertex, 
hep-th/0305132, Commun. Math. Phys. 254 (2005) 425-478. 

[2] A. Aptekarev, P. Bleher, and A. Kuijlaars, Large n limit of Gaussian ran¬ 
dom matrices with external source, part II, math-ph/0408041 

[3] E. Brezin and S. Hikami, Universal singularity at the closure of a gap in a 
random matrix theory, cond-mat/9804023 Phys. Rev. E (3) 57 (1998), no. 
4, 4140-4149. 

[4] E. Brezin and S. Hikami, Level Spacing of Random Matrices in an External 
Source, cond-mat/9804024, Phys. Rev. E (3) 58 (1998), no. 6, part A, 7176- 
7185. 

[5] R. Cerf and R. Kenyon, The low-temperature expansion of the Wulff crystal 
in the 3D Ising model, 

[6] H. Cohn, R. Kenyon, and J. Propp, A variational principle for domino 
tilings, math.CO/0008220, J. Amer. Math. Soc. 14 (2001), 297-346. 

[7] R. Dijkgraaf, A. Sinkovics, and M. Temurhan, Universal Correlators from 
Geometry, hep-th/0406247 

[8] P. L. Eerrari and H. Spohn, Step fluctuations for a faceted crystal, 
cond-mat/0212456, J. Stat. Phys. 113 (2003), 1-46. 

[9] K. Johansson, Discrete polynuclear growth and determinantal processes, 
math.PR/0206208 

[10] K. Johansson, The Arctic circle boundary and the Airy process, 
math.PR/0306216 

[11] V. Kac, Infinite dimensional Lie algebras, Cambridge University Press. 

[12] R. Kenyon, Local statistics of lattice dimers, Ann. Inst. H. Poincare Probab. 
Statist. 33 (1997), no. 5, 591-618. 

[13] R. Kenyon and A. Okounkov, Limit shapes and complex Burgers eguation, 
in preparation. 

[14] A. Okounkov, Random surfaces enumerating algebraic curves, 
math-ph/0412008 

[15] A. Okounkov and N. Reshetikhin, Correlation function of Schur process 
with application to local geometry of a random 3-dimensional Young dia¬ 
gram math.CO/0107056, J. Amer. Math. Soc. 16 (2003), no. 3, 581-603. 

[16] A. Okounkov, N. Reshetikhin, and C. Vafa, Quantum Calabi-Yau and Clas¬ 
sical Crystals, hep-th/0309208 


44 










[17] M. Praehofer and H. Spohn, Scale Invariance of the PNG Droplet and the 
Airy Process, math.PR/0105240, J. Stat. Phys. 108 1071-1106 (2002). 

[18] N. Saulina and C. Vafa, D-branes as Defects in the Calabi-Yau Crystal, 
hep-th/0404246 

[19] C. Tracy and H. Widom, Differential equations for Dyson process, Com- 
mun. Math. Phys. 252 (2003), 7-41. 

[20] C. Tracy and H. Widom, The Pearcey Process, math.PR/0412005 


45 
















